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ABSTRACT 

Recent discoveries by the Rossi X-Ray Timing Explorer indicate that most of the rapidly 
accreting (M ;> 10 _11 M Q yr _1 ) weakly magnetic (B <C 10 11 G) neutron stars in the Galaxy 
are rotating at spin frequencies v s ;> 250 Hz. Remarkably, they all rotate in a narrow range of 
frequencies (no more than a factor of two, with many within 20% of 300 Hz). I suggest that 
these stars rotate fast enough so that, on average, the angular momentum added by accretion 
is lost to gravitational radiation. The strong v s dependence of the angular momentum loss rate 
from gravitational radiation then provides a natural reason for similar spin frequencies. Provided 
that the interior temperature has a large scale asymmetry misaligned from the spin axis, then 
the temperature sensitive electron captures in the deep crust can provide the quadrupole needed 
(~ 10 -7 MR 2 ) to reach this limiting situation at v s rj 300 Hz. This quadrupole is only present 
during accretion and makes it difficult to form radio pulsars with v s > (600 — 800) Hz by accreting 
at M ;> fO _lo M yr _1 . The gravity wave strength is h c ~ (0.5 — 1) X 10~ 26 from many of these 
neutron stars and > 2 x 10~ 26 for Sco X-f. Prior knowledge of the position, spin frequency 
and orbital periods will allow for deep searches for these periodic signals with gravitational wave 
interferometers (LIGO, VIRGO and the "dual-recycled" GEO 600 detector) and experimenters 
need to take such sources into account. Sco X-l will most likely be detected first. 
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1. Introduction 

The launch of the Rossi X-Ray Timing Explorer (RXTE) has allowed for the discovery of fast quasi- 
periodic variability from many rapidly accreting (M ;> 1O -11 M0 yr _1 ) neutron stars. These observations 
strongly suggest that these neutron stars (NSs) are rapidly rotating, as predicted by those scenarios con- 
necting the millisecond radio pulsars to this accreting population (see Bhattacharya 1995 for an overview). 
Strohmayer et al. (1996) were the first to detect nearly coherent vb — 363 Hz oscillations during type I X-ray 
bursts from the low accretion rate (M < 10 _9 M Q yr -1 ) NS 4U 1728-34. Pulsations were detected in six of 
the eight bursts analyzed at that time. In addition, two high frequency quasi-periodic oscillations (QPOs) 
were seen in the persistent emission. These changed with accretion rate, but maintained a fixed difference 
frequency of Vd ~ 363 Hz, identical to the period seen during the bursts. The detection of two drifting QPO's 
(in the persistent emission) separated by a fixed frequency identical to that seen in the bursts naturally leads 
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to beat frequency models (Strohmayer et al. 1996; Miller, Lamb, & Psaltis 1998). The difference frequency 
is presumed to be the NS spin frequency, v s , whereas the upper frequency has different origins in different 
models (see van der Klis 1998 for a summary). In addition, the temporal behavior of the periodic oscillations 
both during the rise of the bursts (Strohmayer, Zhang, & Swank 1997b) and in the cooling tails (Strohmayer 
et al. 1997a) are most easily explained in terms of rotation. 

There are six NSs with measured periodicities during Type I X-ray bursts (see Table 1). Both the 
difference frequencies (vj) and the burst frequencies {vb) are in a narrow range, from 260 to 589 Hz. For two 
objects (KS 1731-260 and 4U 1636-53) the difference frequencies are one-half the burst values. Which value 
is v s is not resolved. There are also many NSs that accrete at higher rates and are not regular Type I X-ray 
bursters. Many of these objects, notably the "Z" sources, also show drifting QPO's at fixed separation, again 
with a similarly narrow frequency range (roughly 250-350 Hz). Beat-Frequency like models are also applied 
to these observations so as to infer v s . The applicability of such a model is less clear when the difference 
frequency is not constant (Sco X-l, van der Klis et al. 1997; 4U 1608-52, Mendez et al. 1998). 

If accreting matter always arrives with the specific angular momentum of a particle orbiting at the NS 
radius (R = 10i? 6 km), then it only takes ~ 10 7 yrs of accretion at M w 10~ 9 M Q yr" 1 for a M = 1AM 1A M Q 
NS to reach v s = 50 Hz from an initially low frequency. It is thus remarkable that these NSs are all rotating 
at nearly the same rate. White and Zhang (1997) argued that this similarity arises because these NSs are 
magnetic and have reached an equilibrium where the magnetosphcric radius equals the co-rotation radius. 
This requires an intrinsic relation between their magnetic dipoles, /if,, and M so that they all reach the same 
rotational equilibrium (most likely \ib oc M 1 / 2 ) and a way of hiding the persistent pulse typically seen from 
a magnetic accretor. 

My alternative explanation for these spin similarities is that gravitational wave (GW) emission has 
started to play an important dnrolc. If the NS has a misaligned quadrupole moment, Q, then the strong 
spin frequency dependence of GW emission defines a critical frequency beyond which accretion can no longer 
spin-up the star. Such a NS will radiate energy via GWs at the rate E = 32GQ 2 o> 6 /5c 5 , where 10 — 2-kv Si and 
lose angular momentum at the rate N gw = E/ui. Balancing this spin-down torque with the characteristic 
spin-up torque from time-averaged accretion, N a w M(GMR) 1 ^ 2 , gives the Q needed so as to make the 
critical frequency 300 Hz, 



or < 10 7 of the NS moment of inertia, / « 10 45 g cm 2 . The similarities in v s may then arise because of the 



What is the source of the misaligned quadrupole? Wagoner (1984) argued that accreting NSs would get 
hung-up at spin frequencies where the Chandrasekhar-Friedman-Schutz (CFS) instability sets in. However, 
Lindblom (1995) and Lindblom & Mendell (1995) have shown that the star needs to be very near the breakup 
frequency (v s ^ kHz) for such an instability to occur, even for the core temperatures T c = (1 — 3) x 10 8 K 
of rapidly accreting NSs (Ayasli & Joss 1978; Brown & Bildsten 1998). The spin frequencies for these NSs 
are too slow for such an instability. 

I present in §2 a new source for lateral density variations in an accreting NS; electron captures (hereafter 
EC) in the crust. The constant compression of the crust forces nuclei to undergo EC when the electron Fermi 
energy, Ef, is high enough to make a transition. However, the crust is hot enough in a rapidly accreting 
NSs to make the EC rates temperature sensitive. Hotter regions then capture at lower pressures, so that 




(1) 



weak dependencies of the critical frequency on Q and M. 
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the density jump associated with the EC is at a higher altitude in the hotter parts of the crust. Moderate 
lateral temperature variations then lead to density variations large enough to generate the required Q. This 
outcome is independent of the particular source of the temperature variations. One possible cause for a T 
asymmetry relative to the spin axis is a weak magnetic field. I conclude in §3 by finding the GW signal 
strength and estimating detection. 



2. Electron Captures in the Neutron Star Crust 

These NSs accrete hydrogen and helium rich material from their companions, and within days of reaching 
the surface, this matter is burned to heavy elements. The composition of these ashes is still uncertain, but 
most certainly consists of heavy nuclei, potentially beyond the iron group (Schatz et al. 1997). This material 
replaces the primordial NS crust and becomes neutron rich by successive EC. Later on, neutron emissions 
and pycnonuclear reactions occur (Bisnovatyi-Kogan & Chechetkin 1979; Sato 1979; Haensel & Zdunik 1990 
(hereafter HZ)). For the M's and residual hydrogen abundance (Taam, Woosley, & Lamb 1996) appropriate 
for these NSs, the temperature of the accumulating matter and deep crust is Tg = T/10 8 K ss (2 — 6) (Brown 
& Bildsten 1998). 

I confine my discussion in this Letter to the outer crust (before neutron drip at p w (4 — 6) x 
10 11 g cm -3 , HZ), which is held up by degenerate and relativistic electrons. They exert a pressure p = 
1.42 x 10 30 erg cm -3 (i?F/30McV) 4 from which the mass of the shell above a particular depth is M cr (E F ) = 
4:irR 2 p/g, where g = GM/R 2 , or 

^)„5xlO-»J»f ^-(^) 4 . (2) 

The crust is compressed on a timescale t comp = p/rhg, where m is the local accretion rate, so that Ep 
rises with time, eventually leading to nuclear EC. Many (Sato 1979; HZ) presumed that the captures were 
instantaneous once Ep > Ed, where Ed is the mass difference between the (A,Z) nucleus and (A,Z-1) nucleus. 
In reality, the transitions will not occur until the EC lifetime is comparable to t cornp (Blaes et al. 1990). For 
low accretion rates M ~ 10~ 16 M Q y _1 , Blaes et al. (1990) showed that most captures occur in a thin zone 
where Ep > Ed only slightly. The situation is different for the high M's of the bright X-ray sources, as the 
high crustal temperature leads to most EC occurring out on the thermal tail at a physical location where 
Ep < Ed- 

I consider a region of the crust consisting of a single nucleus of charge Ze and mass Am p . For many 
nuclei, the capture to the ground state of the (A,Z-1) nucleus is highly forbidden and proceeds more slowly 
than icomp- The element is still abundant when Ep reaches a value where a more favorable transition to an 
excited state of the (A,Z-1) nucleus can occur .[] For this reason, I estimate the EC rates by using ft values 
in the range of allowed reactions, ft = 10 4 — 10 6 s. Since these stars are hot enough for the sub-threshold EC 
to dominate, the transition occurs when Epjepi = Ed — PhsT, where (3 ~ (8 — 12) depends logarithmically 
on m, T and ft (Bildsten & Cumming 1998). A hotter region undergoes the transition at a higher altitude 



1 Since the mass difference of the (A,Z-1) and (A,Z-2) nuclei is smaller than the mass difference of the (A,Z) and (A,Z-1) 
nuclei, the sequence always consists of two successive captures (HZ). This also blocks the inverse reactions. The larger phase 
space for the second reaction (as well as the larger number of available excited states with favorable spin-parity) means that 
the (A,Z-1) nucleus captures electrons more rapidly than the (A,Z) nucleus. Since the first capture is the rate-limiting step, I 
do not track the intermediate (A,Z-1) nucleus. 
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(call s c the distance down to the transition from where Ep — 1 MeV, i.e. s increases into the star) compared 
to a colder region. Figure [l] displays the EC transition layer for a NS accreting at M = 2 x 10 -9 M Q y" 1 . 
As an example I chose the nucleus 56 Ca, which captures at E d ss 23 MeV, where p w 2.5 x 10 11 g cm -3 . 
This is the fourth EC transition that HZ found when starting with 56 Fe. Since this is a generic process, 
I show cases for ft = 10 4 s (dotted) and ft = 10 6 s (solid) for (from left to right in the upper and lower 
panels) Tg — 6,4 and 2. The density jump associated with the ECs is evident and, if discontinuous, would 
be Ap/p = 2/(Z - 2) = 0.11. 
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Fig. 1. — The density, pressure and nuclear abundance in the 56 Ca electron capture layer for a R = 10 km, 
M = 1.4M Q NS accreting at M = 2 x 10 _9 M Q y^ 1 . These are plotted as a function of increasing depth 
into the star, deeper regions are to the right. For a fixed value of ft, the hotter crusts deplete sooner. The 
curves are, from left to right, for T% = 6, 4 and 2. 

If a transverse temperature gradient is present in the crust, then these ECs generate density asymmetries 
on the star and a quadrupole with mass AM « 47ri? 2 ApAs c . The mass inferred from the ft = 10 4 s case 
for a temperature contrast AT§ = 2 (where the difference in physical depth of the EC is As c « 500 cm) 
is AM w 10~ 7 M©, just what is needed for GW's to be important. More generally, the mass involved is 
AM = M cr (Ap/p)(4:AEF idepl /E F ), or since AE F . depl » 10k B AT, 

AM « 5.5 x lO- 8 M AT 8 -^_ ( Ep \ (3) 
° Mi.4 V30 MeV / ' y ' 

for a single EC layer. Since there are a few more electron capture layers underneath this one, there is 
adequate mass to generate the required Q when AT§ « 1. The repercussions of the deeper reactions is still 
unknown. 

The EC's deposit heat directly into the crust, generating a flux Fd ocp = E^cc-pm, where -Edccp ~ 
10 17 erg g _1 (HZ). Temporal variations in rh lead to local heating at these depths, which takes a few 
years to equilibrate in the radial direction. The heat transfer is fixed by electron conduction at these depths, 
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so that the conductivity is K = Tr 2 kgTn e c 2 /ZEpVe-ph, where v e -ph ~ VZaksT /% is the electron-phonon 
scattering frequency (Yakovlev & Urpin 1980). The time it takes heat to cross a scale height, H = p/pg, is 
then i t h s=s pC p H 2 /K ~ yr where C p ~ 3kg /Am p is the specific heat at constant pressure. 

However, lateral thermal asymmetries will persist as long as the heating time, tu ~ C p Ty / Ed ccp rh ~ 
4 x 10 _3 Tgt comp , is shorter than the time to transport heat around the star at the same depth. That time 
is *th,R ~ pC p R 2 /K « 6400 y(PF/30 MeV), so for M ;> 3 x lQ- u M y- 1 (E F /3O MeV) 3 (P/10 km) 2 , the 
time to transport heat around the crust is longer than the time to locally heat it. All of the objects I am 
discussing are in this regime and so large-scale temperature asymmetries can persist as long as the system 
is being perturbed in some way. Though the origin of lateral T gradients is still unknown, it is difficult 
for conduction to wash them out. For example, large-scale m or compositional variations (due to, say, 
magnetically channeled flow or a "buried field" ) will imprint large scale T asymmetries. 

The pressure as a function of depth in the EC region is shown in the middle panel of Figure 1. Even 
though the density contrasts are large, the pressure contrasts are not. If one part of the star has a slightly 
deeper capture zone, at s c ,2 = s c ,i + As c , then the pressures at a fixed depth s > s Cj i relate via P± w 
^2(1 + 8As c /Zs). The pressure is higher underneath the layer that captured sooner due to the "extra 
weight" of the dense layer. Since s cx Ep, the lateral pressure contrast at a fixed depth, AP = (Pi — P2), 
is AP/P w 80/cbAT /ZEp. These transverse pressure gradients would lead to flow and a cancellation of 
the quadrupole if the matter was in a liquid state. However, as long as Tg < 10, the matter is solid at 
these depths and has a finite shear modulus, fi ~ 10~ 2 p (Strohmayer et al. 1991). The transverse pressure 
gradients are then balanced by a slight amount of shear stress as the matter underneath the cooler regions 
laterally shifts a distance £. Roughly, £ is found from the transverse momentum equation AP/P ks p£,/s 2 , or 
£/s ~ (100s/P)(8As c /Zs). This gives (/sw6x 10~ 3 ATg, or a few meters of transverse displacement over 
the whole surface. For now, I will presume that this configuration will not crack. Relative vertical motion 
has yet to be investigated and could reduce the estimated Q. 



3. Detectability of the Gravitational Waves 

Most accreting, weakly magnetic neutron stars are spinning at v a m (250 — 500) Hz. I conjecture 
that this similarity is the result of an equilibrium where the angular momentum accreted is radiated away 
in gravitational waves. I showed that pre-threshold electron captures will turn any lateral temperature 
gradients into lateral density gradients. In this case, a small spin-misaligned temperature gradient gives rise 
to a quadrupole, Q ~ 10~ 7 AfP 2 , adequate to explain the similarities at w 300 Hz. These accretion-induced 
Q's clearly make it difficult to spin up a neutron star to 1000 Hz by accreting at M ;> 10~ 10 A/q yr -1 . Indeed, 
a NS spinning at 1000 Hz with Q w 10~ 8 A/P 2 would be difficult to spin-up even at M ~ 10 -8 Af© yr -1 . 
The lateral T asymmetries should subside at lower Af's, reducing Q (see the dependence in equation ||) 
and potentially making it easier to reach higher v a by accreting slowly for a long period of time. There is 
already some hint of this, as the highest M sources seem to be rotating the slowest. 

Independent of the cause of the T asymmetry, if gravitational radiation is the explanation of the spin 
similarities, then the characteristic GW signal strength at Earth (suitably averaged over spin orientations, 
Brady et al. 1998) is h c = 2.9Gcu 2 Q /dc 4 , where d is the distance to the object (the prefactor is 4 when 
looking down the spin axis). Presuming the NS luminosity is L m GMM/R then h c is written in terms of 
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the observable F = L/And 2 (Wagoner 1984) 

h c w 4 x ltr 27 -^ = ^-r 

M \1" er S cm s 

which represents a lower limit for ft, c due to the L to M conversion I chose and since F is never fully measured. 
The minimum /i c 's for those NSs which have iA, inferred from Type I bursts are shown in Table 1. There 
are many NSs for which only v d has been measured (see van der Klis 1998 for a summary). The average 
2-10 keV fluxes from van Paradijs (1995) imply the following GW amplitudes: GX 349+2 {v d = 266 ± 13, 
h c ~ 5.4 x 1CT 27 , Zhang et al. 1998b), 4U 1820-30 (y d = 275 ±8, h c w 3.7 x 10~ 27 ), GX 17+2 [y d = 294 + 8, 
h c w 4.7 x 10~ 27 ), 4U 0614+09 {v d = 327 ±4, h c » 1.3 x 10~ 27 ), GX 5-1 = 327 ± 11, ft, c « 6.0 x 10~ 27 ) 
Cyg X-2 (y d = 343 ± 21, h c » 3.7 x 10~ 27 ), and GX 340+0 [y d = 325 ± 10, h c » 3.7 x 10^ 27 , Jonker et al. 
1998). 

The highest flux object that shows kHz QPO's is Sco X-l with an orbital period of 18.9 hours. It's 
difference frequency is not constant however (van der Klis et al. 1997) and so we are still uncertain about 
the exact spin frequency. If it is rotating at ~ 250 Hz, then h c « 2.2 x 10~ 26 at 2v s = 500 Hz. Detection 
of a weak periodic gravitational wave signal depends on the ability to coherently fold large data sets for 
a long time, t, in order to obtain the ~ {2v s t) 1 ^ 2 signal enhancement. The initial LIGO interferometer 
will have h rms = 10~ 21 at v « 500 Hz (Abramovici et al. 1992), requiring integration times ~ yr. If 
the orbital parameters, spin frequency and phase were known for Sco X-l, then Brady et al.'s (1998) work 
implies it would be detected at 99% confidence with the initial LIGO in about one year. However, our lack 
of detailed prior information on these parameters demands a larger parameter space search. The consequent 
reduction in the S/N (maybe a factor of ~ 5) make it more likely that detection will have to wait for the 
enhanced LIGO (see Thorne 1998 for the sensitivities). The similar spin frequencies of so many objects 
give an advantage to dual-recycling interferometers (such as GEO 600 and TAMA 300). Indeed, the GEO 
600 noise strength (from Strain et al. 1998) when tuned to Sco X-l is competitive with the initial LIGO 
and VIRGO. Detection of the larger number of NSsQ with h = (0.5 — 1) x 10~ 26 will have to wait for 
advanced broad-band interferometers (such as advanced LIGO), further developments in the dual-recycled 
interferometers, or larger cryogenic detectors. 

The stability of the GW signal depends on the origin of the mass quadrupole. For the EC origin, I would 
not expect major differences in Q on timescales shorter than the thermal time (years). Another relevant time 
is that to replenish the thickness of the EC layer, which is ~ 100 years at M ~ 10 _9 M Q yr -1 . However, since 
M varies, we cannot expect the torques to always be in balance. For example, if M were to drop suddenly, the 
instantaneous spin-down of the NS due to the equilibrium Q would be v w 10~ 13 s _2 ((M)/10 _9 Mq yr -1 ). 
This would lead to one pulse cycle of drift in a time w {2/i>) 1 / 2 w 50 days. 

How does this work affect the accretion spin-up scenario for making millisecond pulsars? The EC Q will 
decrease once accretion has halted and the temperature equilibrates. In addition, when the crust is cold, it 
only takes 10 _7 M Q of accretion to wash out the compositional gradients generated during accretion at high 
rates. There is an important comparison to make in this regard. Namely, as long as Q decays on a timescale, 



) m 



2 The recent RXTE discovery of a persistent 401 Hz accreting pulsar (Wijnands & van der Klis 1998) in the transient binary 
SAX J1808.4-3658 (in't Zand et al. 1998) provides hope for GW studies. It was observed in outburst by RXTE (Marshall 1998) 
and it's orbit measured by Chakrabarty & Morgan (1998). Enough is known about this weakly magnetic (B < 10 9 G) accreting 
pulsar that a coherent data fold can be undertaken. Unfortunately, the time averaged flux (I estimate 2 X 10 — 10 ergs cm -2 s — 1 
at most) is low and even if looking right down the spin-axis, h < 7 X 10~ 28 . It might be that spin-down in quiescence is the 
best indicator of GW emission from this NS. 
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td, much shorter than i S pin-up = Iuji/M(GMR) 1 ^ 2 w 5 x 10 7 yr then the millisecond pulsar is "born" at 
roughly the same spin period since Vf = ^(1 + 4idAspin-up)~ 1 ^ 4 - If Q persists at some low level long after 
accretion has ended, it could still play an important role. Indeed, as many have shown (Thorne 1987, New 
ct al. 1995, Brady et al. 1998) a small Q can go a long way towards explaining the observed spin-down of 
millisecond pulsars and making them detectable as GW sources. 

Transverse T gradients may be more prevalent on the B J> 10 11 G accreting pulsars. However, most 
of them have v s <C Hz so that the Q needed for GW equilibrium is implausibly large, making it unlikely 
that GW's play an important role in their spin evolution. However, there is a lack of accreting pulsars with 
v s > Hz that is usually explained by stating that all pulsars have reached their magnetic equilibrium and 
that few have B <C 10 11 G. If we presume a maximum allowed quadrupole Q max = 10~ 5 MR 2 (Thorne 1987), 
then GW emission could play an important role when v s > 35 Hz(M/10 -9 M Q yr -1 ) 1 / 5 , only somewhat 
alleviating this issue. 

I thank E. Brown, A. Cumming, E. Flanagan, T. Prince, H. Schatz, G. Ushomirsky and I. Wasserman 
for helpful discussions and J. Arons, R. Rutlcdgc, and M. van der Klis for comments on the manuscript. P. 
Brady and K. Thorne clarified issues of GW detection. This research was supported by NASA via grants 
NAG 5-2819 and NAGW-4517, a Hellman Family Faculty Fund Award (UCB) and the Alfred P. Sloan 
Foundation. 
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Table 1: Rapid Periodicities During Type I X-Ray Bursts 



Object Name 






Flux a 


h c 


Ref. b 




(Hz) 


(Hz) 




(io- 27 ) 




4U 1702-429 


330 




1.0 


1.2 


[1] 


4U 1728-34 


363 


363 


2.8 


2.0 


[2] 


KS 1731-260 


524 


260 ± 10 


0.2-2 


2.0 


[3-4] 


Aql X-l 


549 








[5] 


4U 1636-53 


581 


276 ± 10 


4.4 


2.8 


[6-7] 


MXB 1743-29 


589 








[8] 



a Average 2-10 keV fluxes (in units 10 9 erg cm 2 s : ) are from van Paradijs (1995). 

6 REFERENCES: [1] Swank et al. 1997, [2] Strohmayer et al. 1996, [3] Smith et al. 1997, [4] Wijnands & 
van dcr Klis 1997, [5] Zhang et al. 1998a, [6] Strohmayer et al. 1998, [7] Wijnands et al. 1997, [8] Strohmayer 
et al. 1997a 



